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ABSTRACT 

Solar models using the new lower abundances of Asplund et al. (2005, 2009) or Caffau et al. (2008, 
2009) do not agree as well with helioseismic inferences as models that use the higher Grevesse & Noels 
(1993) or Grevesse & Sauval (1998) abundances. Adopting the new abundances leads to models with 
sound speed discrepancies of up to 1 .4% below the base of the convection zone (compared to discrepancies 
of less than 0.4% with the old abundances), a convection zone that is too shallow, and a convection zone 
helium abundance that is too low. Here we review briefly recent attempts to restore agreement, and 
we evaluate three changes to the models: early mass loss, accretion of low-Z material, and convective 
overshoot. One goal of these attempts is to explore models that could preserve the structure in the interior 
obtained with the old abundances while accommodating the new abundances at the surface. Although 
the mass-losing and accretion models show some improvement in agreement with seismic constraints, a 
satisfactory resolution to the solar abundance problem remains to be found. In addition, we perform a 
preliminary analysis of models with the Caffau et al. (2008, 2009) abundances that shows that the sound 
speed discrepancy is reduced to only about 0.6% at the convection zone base, compared to 1.4% for the 
Asplund et al. (2005) abundances and 0.4% for the Grevesse & Noels (1993) abundances. Furthermore, 
including mass loss in models with the Caffau et al. (2008, 2009) abundances may improve sound speed 
agreement and help resolve the solar lithium problem. 

Subject headings: Sun: abundances - Sun: evolution - Sun: helioseismology - Sun: interior - Sun: oscillations 



1. Introduction 

1.1. The Solar Abundance Problem 

Before 2005, we believed we knew very well how 
to model the evolution and interior structure of the sun. 
Evolved solar models with the latest input physics (in- 
cluding diffusive helium and element settling and with- 
out tachocline mixing) reproduced the sound speed 
profile determined from seismic inversions to within 
0.4%, as well as the seismically-inferred convection 
zone depth and convection zone helium abundance. 
Solar interior modelers had little impetus to progress 
beyond one-dimensional spherical models of the sun, 
with perhaps the exception of introducing some ad- 
ditional mixing below the convection zone (hereafter 
CZ) to deplete surface lithium and reduce the small 
remaining sound-speed dis c repanc y at the CZ base. 
For the brevesse & Noelsl (1 19931 hereafter GN93) 



or iGrevesse & Sauvall d!998l hereafter GS98) abun- 
dances, the simplest 'spherical sun' assumptions ap- 
peared nearly adequate for solar modeling. These in- 
clude one-dimensional zoning (concentric shells in 
hydrostatic equilibrium), initial homogeneous compo- 
sition, negligible mass loss or accretion, neglecting 
rotation and magnetic fields, simple surface bound- 
ary c onditions, mix i ng-len gth theory of convection 
(e.g., iBohm-Vitensd (119581) ). and no additional mix- 
ing or structural changes from convective overshoot, 
shear from differential rotation, meridional circulation, 
waves, or oscillations. 

However, new analyses of solar spectral lines re- 
vise downward the abundances of elements heavier 
than hydrogen and helium, particularly the abundances 
of carbon, nitrogen, and oxygen that contribute to the 
opacity just below the CZ. See Table Q]for a summary 
of some of the major abundance revisions over the last 
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Table 1 : Mass fractions of metals in the present-day photosphere, Z, and ratio of metals to hydrogen mass fraction, 
Z/X, evaluated over the last twenty years. 



Year 


Source 


Z 


Z/X 


1989 


Anders & Grevesse (1989) 


0.0201 


0.0274 


1993 


Grevesse & Noels (1993) 


0.0179 


0.0244 


1998 


Grevesse & Sauval (1998) 


0.0170 


0.0231 


2005 


Asplund et al. (2005) 


0.0122 


0.0165 


2009 


Asplund et al. (2009) 


0.0134 


0.0181 


2009 


Ludwig et al. (2009) 


0.0154 


0.0209 



"Values for Z are inferred from Z/X assuming Y=0.248. 
''Uncertainties are <10%. 



twenty years. 

Compare d to the older GS98 and GN93 abun- 
dances, the lAsplund etail d2005[ hereafter AGS05) 
abundance of C is lower by 35%, N by 27.5%, O by 
48%, and Ne by 74%. The abundances of elements 
from Na to Ca are lower by 12 to 25%, and Fe is de- 
creased by 12%. For the GS98 abundances, the ratio 
of the element to hydrogen mass fraction Z/X = 0.023, 
and the heavy element abundance Z ~ 0.018, while, 
for the new abundances, Z/X = 0.0165, and Z~0.0122. 
Models evolved with the AGS05 abundance mixture 
give worse agreement with helioseismic constraints; 
the sound-speed discrepancy is 1.4% below the CZ 
base, and the CZ depth is shallow and CZ helium abun- 
dance is low compared to those derived from seismic 
inversions. 

Recently. lAsplund et alj d2009l hereafter AGSS09) 
re-evaluated the spectroscopic abundances, carefully 
considering the atomic input data and selection of 
spectral lines and using improved radiative trans- 
fer and opacities. They revised the heavy element 
abundance to Z/X = 0.0181 and Z = 0.0134. This 
slight increase over the AGS05 values yields some 
impr ovement in agreeme nt with seismic constraints 
(see Serenelli et al. I l2009h . although the highei abun- 
dances of GN93 or GS98 still provide the best agree- 
ment. In addition, the Cosmological Impact of the 
First STars Team and its collaborators used the 3D 
model atmosphere code CO 5 BOLD to perform an 
independent investigation of the solar abundances 



dCaffau et al.ll2008albl 120091: iLudwig et al.ll2009l here- 
after CO s BOLD). They derived heavy element abun- 
dances of Z/X = 0.0209 and Z = 0.0154, in between 
the AGSS09 and GN93 values. 



Spectroscopic determinations measure the pho- 
tospheric abundances. The continuous convective 
overshoot into the photosphere should leave the 
photosphere with the same abundances as the con- 
vection zone. In addition, convective timescales 
are much shorter than element diffusion and evo- 
lution timescales, making the convection zone well 
mixed and homogeneous. Therefore the spectroscopic 
abundances should be indicative of the abundances 
throughout the convection zone. However, we are re- 
luctant to dismiss the recent abundance re-analyses 
because of the many improvements in the physics and 
models included, namely 3D dynamical atmosphere 
models, non-local thermodynamic equilibrium correc- 
tions for important elements, and updated atomic and 
molecular data. Line profile shapes now agree nearly 
perfectly with observations. Also, it is impressive that 
abundances derived from several different atomic and 
molecular lines for the same element now are consis- 
tent. 

For this paper, we will focus on the AGS05 abun- 
dances as the most extreme example of the lower abun- 
dance determinations. However, we will devote sec- 
tion [7] to a preliminary exploration of the CO 5 BOLD 
abundances. We first review the results of helioseis- 
mic tests using the old and new abundances and on- 
going attempts to resolve these discrepancies (Section 
Q]). We discuss in detail the following three mitigation 
attempts: an early mass-loss phase in solar evolution 
(Section [3]), accretion of low-Z material early in the 
sun's lifetime (Section [H, and extending the CZ be- 
low the depth inferred by helioseismology (Section|5]). 
We then examine including some of these adjustments 
in models with the higher metallicity of CO 5 BOLD 
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(Section |7). The models used to explore all of these 
changes are described in Section[2] 

1.2. Helioseismology and solar models 

Helio- and asteroseismology have turned out to be 
excellent tool s to te s t the physics of stellar models . 
Basu & Antial d2004;lBahcall & Pinsonneaultl d2004 : 
Turck-Chieze et al.1 (120041) authored some of the first 
papers to examine the effects of the new lower abun- 
dances on solar models, demonstrating that the new 
abundances lead to greater discrepancy with seismic 
inferences. This is demonstrated in Table [2] which 
compares our calibrated evolution models using the 
GN93 and AGS05 mixtures (details of our models can 
be found in Section|2]i. For the AGS05 model, the CZ 
helium abundance Y is low (0.22730) and the CZ base 
(Rczb = 0.72944 R ) is shallow compared to the seis- 
mically inferred CZ Y abundance of 0. 248±0.003 and 



CZ base radius of 0.713+0.001 R n from Basu & Antia 



(12004 . 

Figure[TJshows the differences between inferred and 
calculated sound speed for calibrated evolved mod- 
els using the old and new abundances. The inferred 
sound speed is from iBasu et al. ( 2000h . The uncer- 



tainties in sound speed inversions are much smaller 
than the differences between these curves, at most a 
few widths of the plotting line. Figure [2] shows the 
observed minus calculated frequency differences for 
modes of angular degree t = 0, 2, 10, and 20 that 
propagate into the solar interior below the convection 
zo ne. The calcula ted frequencies were computed using 
the Pesnell (1990) non-adiabatic stellar p ulsation code 
The o bservatio nal data are from BiSON (jChaplin et al. 
20071) LowL dSchou & Tomczykl fl996l) . or GOLF 



dGarcia et alJl200llK The observational uncertainties 
for the modes are less than 0. 1 juHz, much smaller than 
the mod el discrepancies . The O -C trends for the model 
with the lFerguson et al. (2005) low-temperature opac- 
ities are flatter for higher frequencies that are more 
sensitive to the solar surface. These newer opacities 
include an improved treatment of grains, finer wave- 
length spacing, and additional molecular lines, and are 
higher by 12% up to a factor of three. As illustrated by 
both of these plots, the discrepancy with the new abun- 
dances is much larger than with the old abundances. 

Helioseismology has als o been used to inve stigate 
heavy element abundances. iLin & Dappenl d2005l) find 
that a reduction in carbon abundance, in the direc- 
tion of the AGS05 abundances, can improve sound 
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Fig. 1. — Difference between inferred and calculated sound 
speeds with error bars for models with the GN93 and AGS0 5 abun - 
dances. The sound speed inversion is from IBasu et all (12000). 
The seismically infe rred convection zone base at R = 0.713 Rq 
(Basu & Antidl2004h is shown with the vertical line. 
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Fig. 2. — Observed minus calculated versus calculated frequency 
for our models with the GN93 mixture (triangles or circles), and 
the AGS05 mixture (diamonds) for modes of degree t = 0, 2, 10, 
and 20. T he models producing th e triangle and diamond points use 
the newer Ferguson et al. (2005) low-temperature opacities, while 
the model with the circle points uses 1995 Alexander & Fergu- 
son opacities. The calculated frequencies were computed using the 
iPesnel l 11990) non- adiaba ti c stellar pulsation code, and the data are 
fromlChaplin et all fe007l) : ISchou & Tomczykl <1996l) ; lGarcia et all 
(2001). 
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speed inference. However. iLin et al.1 (12007) show that 
lower Z increas es the discrep a ncy w ith adiabatic in- 
dex inversions; lAntia & Basul (120061) use the ioniza- 
tion signature in the sound speed derivative to infer 
Z cz = 0.0172+0.002. These results favor the old, 



of the CZ base. The structure of most of the con- 
vecti on zone is essentially independent of the opac- 



higher abundances. In addition. IChaplin et al.l (2007) 
use small frequency separations between low-degree 
modes that are sensitive to the core structure to con- 
strain the core Z abundance. They find Z core = 0.0187- 
0.0239. IZaatri et alj d2007l) find that the mean low- 
degree frequency spacings of a model using AGS05 
abundances are incompatible wi th those dete r mined 



from the GOLF measu rements of Lazrek et al. (2007) 
and iGelly etal] d2002h . SimilarlyT lBasu et all 120071) 
find that models constructed with low metallicity are 
incompatible with the small frequency spacing and fre- 
quency separation ra tios calculated from BiSON data 
(IChaplin et al.lll996l) . These results do not rule out 
accretion, enhanced diffusion, or other options that 
can retain high core Z. However, they do disfavor the 
prospect that the new, lower abundances were present 
initially throughout the sun. 



1.3. Attempts to restore agreement through solar 
model modifications 

Here we briefly review recent attempts to adjust so- 
lar models in order to mitigate the discrepancy with 
seismic constraints for the new abundances. For a 
more d etailed review of previ ous attempts, see, for ex- 
ample, Basu & Antial (|2008l). In the following sec- 
tions, we discuss increased opacities below the CZ 
(11 - 30%); increased neon abundance (x~4); in- 
creased abundances (within uncertainty limits, or us- 
ing alternative determinations); enhanced diffusive set- 
tling rates (xl.5 or more); accretion of lower-Z ma- 
terial early in the sun's lifetime; structure modifica- 
tion below the CZ base due to radiative damping of 
gravity waves; tachocline mixing (also used with old 
abundances); convective overshoot; and combinations 
of the above. The conclusion has been that it is difficult 
to match simultaneously the new Z/X and helioseis- 
mic constraints for CZ depth, sound speed and density 
profiles, and CZ helium abundance by applying these 
changes. 



1.4. Opacity increases 

Heavy-element abundances primarily affect solar 
structure through their effect on opacity, which affects 
the structure of the radiative zone and the location 



ity. IChristensen-Dalsgaard et al.1 (120091) determined 



the change in opacity required to restore the sound- 
speed agreement of a solar model using the AGS05 
abundances to the level of success originally attained 
with the GN93 abundances. They find that opacity 
would need to be increased ranging from about 30% 
below the CZ to a few percent in the core. 

Although improvements to the solar model can 
be made by ad hoc opacity increases, there is lit- 
tle justification for such large enhancements. The 
presently available opacity tables from three sepa- 
rate projects (OPAL, OP, and LANL T-4) for con- 
ditions below the CZ differ by only a few percent 



dNeuforee-Verheecke et al.H200UlBadnell et al.ll2005l) . 
making it difficult to justify such large opacity en- 
hancements. Using the Los Ala mos National Lab- 
oratory T-4 opac i ty library data dMagee et al.l Il995 



Hak el et alJl2006b . lGuzik et al.l d2009h find that, to ob 



tain a 30% opacity increase with the new abundances, 
the contribution from oxygen alone would need to in- 
crease by a factor of two to three. Alternatively, the 
iron absorption contribution would need to increase by 
a factor of three. 

Including additional elements has a negligible ef- 
fect on Rosseland mean opacities for solar interior 
conditions. The Lawrence Livermore OPAL opacities 
for the AGS05 mixture included 17 of the most abun- 
dant elements. W ith the LANL T-4 opacity library, 
Guzik et al.l (120091) find that including in the mixture 
all of the elements up to atomic number Z=30 in- 
creases the mixture opacity by only 0.2% for solar in- 
terior conditions. Including additional elements in the 
mixture from 30 < Z < 93, an 83-element mixture, fur- 
ther increases the mixture opacity by less than 0.1%. 
An as-yet unidentified error in calculations of the line 
wings for K-shell transitions in O, C, N, and Ne at en- 
ergies around 800 eV, the peak of the weighting of 
the Rosseland mean opacity for the temperature con- 
ditions at the CZ base, could provide some increase in 



the ca lculated Rosseland mean opacities dGuzik et al 
20091) . 



In o rder to shed light on the issue. lTurck-Chieze et al 
(120091) suggest experimental investigations of opac- 
ity coefficients for the radiative zones of solar-like 
stars. The coming large laser facilities (LIL +PETAL, 
OMEGA EP, FIREX II, LMJ, NIF) have the potential 
to attain sufficiently high temperatures and high densi- 
ties at LTE along with the precise diagnostics required 
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for stellar opacity measurements. Experiments to in 
vestigate the opacities relevant for stellar interiors are 



being conducted at Sandia's Z facility by lBailey et al, 
(2009). 



1.5. Neon and other element abundance increases 

For a while, it was thought that an increase in the 
solar neon abundance provided the most plausible res- 
olution to this problem. Neon is not measured in the 
photosphere due to a lack of suitable spectral lines. 
Instead, its abundance is determined relative to oxy- 
gen using lines formed in the solar corona, XUV and 
gamma ray spectroscopy of quiet and active regions, 
and solar wind particle collections. AGS05 adopt a 
Ne/O abundance ratio of 0.15, and apply this ratio 
to the photospheric oxygen abundance to derive the 
neon abundance. The neon abundance has been re- 
vised downward by 74% from the GS98 value, for the 
most part due to the oxygen abundance reduction. 

Several groups explored solar models with en- 
hanced neon. Some improvement in agreement was 



( Antia & Basul2005:iBahcal 


et al.l2005al Turck-Chieze et al. 


20051 IZaatrietal. 2007: 


Delahave & Pinsonneaultl 


2006). However. iLin et al. (|20O7|) find that increas- 



ing Ne alone actually increases the discrepancy in the 
adiabatic exponent in the region 0.75-0.9 R Q . They 
find that the discrepancy is reduced if only C, N, and 
O abundances are increased. 

More modest Ne enhancements, combined with in- 
creases in the other element abundances of ~0.05 dex, 
at the limit of the AGS05 uncert ainties, have a l so been 
considered. The best model of Bahcall et al. (2005a) 
has Ne enhanced by 0.45 dex (2.8x), Ar by 0.4 dex, 
and C, N, and O by 0.05 dex. This model produces rea- 
sonably good agreement with the inferred sound speed 
and density profile, and has CZ base radius 0.715 R 
and acceptable CZ Y= 0.2439. Of course any increase 
in abundances from the AGS05 value will mitigate the 
problem by increasing opacitie s be low the CZ. For ex 



ample s, see IZaatri et al.1 (120071) or lTurck-Chieze et al 
(200i). 



1.6. Enhanced diffusion 

Sev e ral groups, e. g., lBasu& Antial(l2004 . Montalban et al 
(120041). iGuzik et alJ 12005] hereafter GWC05), and 
lYang & Bil (120071) considered the effects of enhanced 
diffusion. At first this idea might seem promising, be- 
cause the solar interior could have higher abundances 



that give good sound speed agreement, while the CZ 
elements could be depleted to the AGS05 photospheric 
values. In practice, the required diffusion increases are 
quite large (factors of 1 .5 to 2 on absolute rates), and 
enhanced diffusion also depletes the CZ Y abundance 
to well below the seismic determination and leaves the 
CZ too shallow. 

GWC05 investigate the enhancement of thermal 
diffusion for elements and He by different amounts. 
A model with resistance coefficients xl/4 for C, N, O, 
Ne, and Mg and x2/3 for He shows some improvement 
in the sound-speed discrepancy, but the CZ depth is 
still a little shallow (0.718 R ), and the CZ Y is still 
a little low (0.227). Moreover, there is no justification 
for these ad hoc changes in thermal diffusion coeffi- 
cients. The diffusion coefficients themselves should not 
be in error by such a large factor. 

1.7. Gravity waves and dynamical effects 

Arnett, Meakin, & Young (2006, privat e communi- 
catio n) have been inve s tigatin g, following lPress ( 1981 ) 
and iPress & Rvbickil dl98ll) . the effects of gravity 
waves excited and launched inward at the CZ base. 
The radiative damping of these waves as they travel 
inward deposits energy and changes the solar structure 
in the same way as would an opacity enhancement. 
The amount of damping and distance that the waves 
propagate depend on the initial amplitudes and the 
degree of the mode, with low-frequency, high-degree 
waves damped more heavily, after traveling a shorter 
distance. The expected wave spectrum and amplitudes 
still need to be worked out, but could remove as much 
as one third of the sound- speed discrepancy. More re- 
cently, lArnett et al.1 (120091) have re-examined convec- 
tion at the surface and sub-surface layers of the sun, 
proposing a way to eliminate astronomical calibration 
from stellar convection theory. By choosing charac- 
teristic lengths that are determined by the flow, they 
eliminate the need for the free parameters traditionally 
used in mixing length theory. They show that some 
of the discrepancy between the new abundances and 
helioseismic inferences may result from the neglect of 
hydrodynamic processes in the standard solar model. 



1.8. Combinations of effects 

In addition to single changes to solar models, 
several groups considered combinations of changes, 
such as diffusion, opacity, an d abundance enhance - 
ments. See, for example, iBasu & Antial (|2004), 
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Montalban etail d2004 . and lBahcall et all d2006l) . 



Although the above modifications to the input 
physics have achieved some success in restoring agree- 
ment between seismic constraints and models that use 
the new abundances, agreement is not fully restored 
and there is little physical justification for the proposed 
changes. In this paper, we discuss the motivation for 
and the results of three additional attempts to restore 
agreement: an early mass loss phase, accretion of low- 
Z material, and convective overshoot. 

2. Solar evolution models 

Sol ar models require input data for o pacities such as 



ratio a are adjusted so that the final luminosity, radius, 
and surface Z/X match the observation a l con straints to 
within uncertainties. See iGuzik et al.l (12005b for ref- 
erences and a description of the physics used in the 
evolution and pulsation codes and models. 

3. Mass loss 

3.1. Motivation and method 



Willson et al.l (11987) explored the possibility that 



OPAL (Iglesias & Rogers 1996) or OP (Seaton & Badnell 



ipa 
(S. 



2004) supplemented by l ow-temperature opacities 



(e.g., Ferguso n et al.l 2005 ); equation of state such 
as OP AL dRogersetal.ll 1996b. MHD dDappen et al ' 



1988), or CEFF (Christensen-Dalsgaard & Dappen 
1992 ): nuclear reaction rates (e.g. JAngulo et a" 1999b: 



and diffusive element settling (e .g., iBurgers 



Cox et al]ll989tlThoul et all l994). 



1969; 



The solar models produced by the Los Alamos 
group shown hereare evolved from the pre-main se- 
quence using an updated versio n of the one-dimensio nal 
evolution codes described in llbenl (1963, 1965 allbh. 



significant mass loss could occur during the early 
part of the main-sequence for ~ 1-2.5 M stars. They 
considered mass-loss rates ranging from 10~ 9 to 10~ 8 
M /yr that would remove a substantial fraction of mass 
from a star before it evolves off the main sequence. 
Mass loss in these stars, possibly including the early 
sun, would be driven by pulsation, which provides 
the necessary mechanical energy flux, and facilitated 
by rapid rotation. The mass-loss rate would diminish 
upon the development of a surface convection zone, 
which channels mechanical energy away from pulsa- 
tion, and of magnetic fields which provide a ngular mo 



mentum transfer and rotational braking. Guzi k et al 



d 19871) showed that mass-losing solar models have 
steeper molecular weight gradients, shorter main- 
.sequence lifetimes, higher 8 B neutrino fluxes, deeper 
The e volution code uses the SIREFF EOS (see lGuzik & Swensog[ lr f ace convection, higher surface 3 He abundances, 



1 1997b. Burgers' d iffusion treatment as implemented 
bv ICox et al dl98S ), the nuclear reaction rates from 
En gulo et al .1 (1999) with a correction to the 14 N rate 



Formicola et al. (2004), and the OPAL opac- 



from 

ities ( Iglesias & Rogers! 1996) supplemented by the 
IFerguson et al.l d2005b or Alexander & Ferguson (pri- 
vate communication , 1995) low-temperature opacities. 
The IFerguson et al. I d2005h low-temperature opacities 
include an improved treatment of grains, finer wave- 
length spacing, and additional molecular lines and are 
higher than the Alexander & Ferguson (1995) low- 
temperature opacities by 12% up to a factor of three. 
As discussed in Section 11.21 the observed minus cal- 
culated frequency trends for a model with the newer 
low-T opacities are flatter for higher frequencies that 
are sensitive to the solar surface, as seen in Figure|2] 

The models ar e calibrated to the present solar radius 
6.9599X10 10 cm dAllenll 1973b. luminosity 3.846xl0 33 _served and calculated />mode oscillation frequencies 
erg/gdWillson et al.|l986|), mass 1.989xl0 3 3 g dCohen & TavlorL test me structure of solar models including early 
Il986l) . age 4.54+0.04 Gyr dGuenther et alJI 1992b . 

an d main-sequence mass loss. They show that extreme 

adopted photospheric Z/X ratio. For evolution models, ^lar mass \ oss has a significant effect on solar struc- 

the initial helium abundance Y, initial element mass ture and can be ruled out by the p-mode oscillation 

fraction Z, and mixing length to pressure-scale-height frequencies. 



and earlier, more pronounced dredge-up of CN cycle 
processed material in the post-main-sequence phase 
compared to the standard model. In addition, mass- 
losing models predict complete destruction of pro- 
tosolar Li and Be, requiring a mechanism, such as 
production in spallation reactions or flares, for partial 
replenishment to the observed surface abundances. 

Such mass loss in other stars could potentially ex- 
plain blue stragglers and the earlier-than-predicted 
dredge-up of carbon and nitrogen in so lar-mass star s 
ascending the first red giant branch (see lGuziklll988b . 
However, there are also drawbacks. If the sun re- 
mains at too high a mass for too long, all surface Li 
in the subsequent solar model at present age is de- 
stroyed, too much surface 3 He is produced, a nd dis- 
crepa ncies with the inferred sound speed arise ( Guzikl 



1988). In addition. IGuzik & Coxl (119951) compare ob- 
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Minton & Malhotra 




While more extreme early main-sequence mass loss 
has not been observed, smaller mass loss in the sun of 
about 0. 1 M looks promising to solve a number of 
problems. The advantages of an early mass loss phase 
in solar evolution include: solving the faint early sun 
problem, explaining early liquid water on Mars, early 
inner solar system bombardment, and solar lithium de- 
struction. 

Models with early mass loss using the older, 
high er element abunda n ces were explored previousl y 
bv IGuzik et all d 19871) . ISwenson & Faulknerl d 19921) 
iGuzik & Col3 dl995h ._and_ 
( 2003b. In addition, 
cently assessed consequences for Earth climate and 
solar system formation. Here we re-investigate so- 
lar mass loss in light of the new abundances. In our 
m odels, we use th e mass- loss treatment impl emented 
bv lBrunishl dl98ll) in the llbenl d 19631 1 1965 alibi) code. 
We evolve two models with initial masses 1.3 and 
1.15 M , having exponentially-decaying mass-loss 
rates with an e-fol ding time of 0.45 Gyr. Following 
Guzik et alj d 1987b . we adopt this exponential mass- 
loss prescription because it is simple and decreases 
smoothly with time. This is a physically plausible de- 
scription, as the mass-loss rate should be highest when 
a rotating star arrives on the main sequence within the 
instability strip, where pulsation and rotation can fa- 
cilitate mass loss. The mass loss should then decrease 
as the star moves out of the instability strip, ceases 
to pulsate, spins down, and develops a surface con- 
vection zone. The present solar mass -loss rate is 2 
xl0~ 14 M /yr (e.g jFeldman et alj[l977l) . too small to 
affect the sun's evolution. The initial mass-loss rates 
for the two models are 6.55 and 3.38 xl0~ 10 M /yr, 
respectively. 
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3.2. Results 

Table|2]summarizes the initial Y and mixing-length 
parameter needed to calibrate each model and the final 
CZ Y and CZ base radius. The seismically-inferred CZ 
Y abundance and CZ base radius a re 0.248 + 0.003 an d 
0.713 + 0.001 R , respectively dBasu & Antiall2004 . 



Figure [3] shows the luminosity versus time for these 
models, as well as for two constant, one solar-mass 
calibrated models. Figure |4] shows the inferred mi- 
nus calculated sound speed for these models. For the 
models with the AGS05 abundances, the sound speed 
agreement is considerably improved by including early 
mass loss. For the model with initial mass 1.3 M 
sound-speed agreement is almost restored near the CZ 
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Fig. 3. — Luminosity versus time for standard one solar-mass 
models using the GN93 and AGS05 abundances and for two mass- 
losing models using the AGS05 abundances with initial mass 1.3 and 
1.15 M Q . Mass-loss rates are exponentially decaying with e-folding 
time 0.45 Gyr. 
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Fig. 4. — Inferred minus calculated sound speed differences 
for calibrated standard one solar-mass models using the GN93 and 
AGS05 abundances and for models with AGS05 abundances and 
initial mass 1.3 and 1.15 M Q including early mass loss. 



base, but the agreement is not as good in the more H- 
depleted core. Unfortunately, while the model with 
initial mass 1.15 M has a little better sound speed 
agreement in the central 0.1 R , the improvement is 
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cies for calibrated standard one solar mass models using the GN93 
and AGS05 abundances, and for models with AGS05 abundances 
with initial mass 1.3 and 1.15 M including early mass loss. Fre- 
quencies compared are fo r modes of angular degrees I = 0, 2, 10, 
and 20 . The data are from Chaplin et al. ( 2007), Scho u & Tomczvlj 
(1996), and Garcia et al. l 200l [). T he calculated frequencies were 
computed using the Pesnell 1 1990) non-adiabatic stellar pulsation 
code. 



Time (Gyr) 

Fig. 6. — Temperature experienced by the present-day solar sur- 
face layer as a function of time for the mass-losing and standard 
models. For the mass-losing phase, the lithium-destroying temper- 
atures are attained because the layer that is now at the surface once 
resided deeper inside the sun. In the post-mass-loss phase, the rele- 
vant temperatures are attained by envelope convection which mixes 
surface layers downward, exposing the surface material to the tem- 
perature at the CZ base. 2.8 million K is the temperature required 
for relatively rapid Li destruction. 



not as pronounced for the region below the CZ. Fig- 
ure [5] shows the observed minus calculated versus cal- 
culated nonadiabatic frequencies for modes of angular 
degrees I = 0, 2, 10, and 20 that propagate into the so- 
lar interior below the convection zone. Including mass 
loss improves agreement, but models with old abun- 
dances and no mass loss still give the best agreement. 

The mass-losing models described here would de- 
stroy all of the observed surface lithium. Lithium 
is destroyed relatively rapidly in the solar interior at 
temperatures > 2.8 million K. For standard models, 
on the main sequence the surface layers are never 
mixed to high enough temperatures to deplete Li by 
the observe d factor of 150 from the initial solar system 
abundance (|Asplund et al] f2009). and additional mix- 
ing mechanisms must be invoked. However, with mass 
loss, layers that are now at the surface were initially in 
the interior at temperatures high enough to quickly de- 
stroy Li. Figure [6] shows the temperature experienced 
by the surface layer throughout the evolution of each 
model. During the mass-losing phase, the high temper- 
atures that depleted the lithium in the current surface 
layer were experienced when the material that is at the 
surface of the now 1 M sun was deeper, before the 



previous surface layers were lost. After the mass-loss 
phase, the temperature that affects the surface layer is 
that of the CZ base, since the material currently at the 
surface is continually mixed through the CZ. 

The mass-losing models also produce more 3 He 
at the surface as the now-surface layers were once 
processed at higher interior temperatures where 3 He 
builds up to higher equilibrium values. For the 1.3 M 
initial mass model, the surface 3 He mass fraction is en- 
hanced from its initial value of 5.0 xl0~ 5 to 9.0 xlO" 5 , 
while for the 1.15 M initial mass model, the surface 
3 He mass fraction is only slightly enhanced from its 
initial value of 5.0 xl0~ 5 to 5.1 xl0~ 5 . The final 
3 He/ 4 He abundance ratios for the 1 .3 and 1.15 M Q ini- 
tial mass models are 3.9 xl0~ 4 and 2.2 xl0~ 4 , respec- 
tively, enhanced from an initial value of 2.0 xlO -4 . 

4. Accretion of low-Z material 
4.1. Motivation and method 

As a way to keep the solar interi or more like models 
obtained using higher abundances. iGuzik et alJ (12004 , 
hereafter GWC04) and GWC05 proposed accretion of 
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material depleted in heavier elements early in the sun's 
lifetime. In this scenario, the pre-main sequence sun 
would have ~98% of its present mass and a higher Z 
with a mixture similar to the GN93 or GS98 abun- 
dances. After the sun begins core hydrogen burn- 
ing and is no longer fully convective, the remaining 
~2% of material accreted would have lower Z provid- 
ing a convection-zone abundance similar to the cur- 
rent photospheric abundances of AGS05 or AGSS09. 
Po ssible justificatio ns fo r this scenario are dis cussed 
bv lNordlundl d2009l) and lMelendez et al.1 d2009l) . One 
plausible explanation is that planet formation removes 
some high-Z elements from the solar nebula, leaving 
lower-Z material to be accreted after the sun is no 
longer fully convective. 

Winnick et al. d2002b explored the accretion of 



metal-rich material in models with the older GS98 
abundances. They show that some solar models with 
enhanced metallicity in the convection zone might be 
viable as small perturbatio n s to the standard GS98 



model. lHaxton & Serenellil (12008b discuss the ques- 



tion of accretion of metal-depleted gas onto the sun as 
a motivation for future experiments to measure CN- 
cycle neutrinos. The flux of these neutrinos should 
depend nearly linearly on the initial core abundance of 
C and N. A successful measurement of the CN-cycle 
solar neutrino flux would therefore place constraints 
on possible accretion by determining the metallicity of 
the solar core. 

To test the possibility of low -Z accretion in the early 
sun, a model is evolved starting with Z = 0.0197 on 
the zero-age main sequence, and material is progres- 
sively added to reduce the CZ Z by 0.001 in each of 
six steps of about six million years. After each step, 
the model is given time to equilibrate to a new shal- 
lower CZ dept h, leaving behind the higher-Z composi- 
tion gradient ( Guzikl2006l) . The final accretion episode 
leaves the CZ with Z = 0.0137. After 36 million years 
of low-Z accretion, the model is evolved normally (in- 
cluding diffusive settling) and calibrated as usual to the 
observed luminosity, radius, and AGS05 Z/X value. 
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Fig. 7. — Heavy element abundance profiles for the GN93 mixture 
model and the low-Z accretion model. 



> 











GN93 

Accretion 























































































0.296 
0.288 
0.28 
0.272 
0.264 
0.256 
0.248 
0.24 

0.2 0.4 0.6 0.8 1 

R(R ) 

sun 

Fig. 8. — Helium abundance profiles for the GN93 mixture model 
and the low-Z accretion model. 



4.2. Results 

Figure [7] shows the heavy element abundance 
throughout the sun in the accretion model, and Figure 
|8]shows the Y abundance. As intended, the abundance 
in the interior is similar to the GN93 model while the 
CZ abundance matches the new, lower Z from photo- 
spheric observations. 



Figure|9]shows the relative sound-speed differences 
for models with the GN93 mixture, the AGS05 mix- 
ture, and low-Z accretion. The accretion model shows 
improvement in sound-speed agreement in the interior 
where Z is similar to the GN93 mixture. However, 
discrepancy remains near the CZ base. Compared to 
the AGS05 model, the accretion model has a less shal- 
low CZ base radius of 0.7235 R and a nearly accept- 
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Fig. 9. — Relative difference between inferred and calculated 
sound speeds for models with the GN93 and AGS05 abundances 
and with low-Z accretion. 



able CZ Y abundance of 0.2407. Figure [10] shows 
the observed minus calculated versus calculated non- 
adiabatic frequencies for modes of angular degrees ( 
= 0, 2, 10, and 20 that propagate into the solar interior 
below the CZ. Including accretion in a model using 
the new abundances improves agreement with this fre- 
quency data. 



Castro et al.l (12007b also approximated an accretion 
model by instantaneously decreasing the Z abundance 
gradient in the CZ in an early main-sequence model 
(age 74 My). They do not find an improvement in the 
CZ depth, as we did for our model, but find about the 
same CZ Y abundance, 0.240, and improved sound- 
speed agreement below 0.5 R Q . 

5. Convective overshoot 
5.1. Motivation and method 

It is possible that the CZ depth predicted using 
standard mixing-length theory is too shallow and con- 
vective motions extend the nearly adiabatically strat- 
ified par t of the CZ to the depth inferred seismically. 
Rempel (2004) presents a semi-analytical model of 
overshoot. This approach facilitates the understand- 
ing of the relation between numerical simulations and 
classical o vershoot theories i n term s of physical pa- 
rameters. iMontalban et al. ( 20061) developed solar 
models that include convective overshoot. By adopt- 
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Fig. 10. — Observed minus calculated frequency versus calcu- 
lated ferquency of GN93, AGS05, and accretion models for de- 
gree 1=0, 2, 10, an d 20 modes. T he calculated frequencies were 
computed using the lPesnell U99QD non-adiab at ic stellar pulsation 
code. The data are from [Chaplin etal J <2007l) , ISchou & Tomczv^ 
<199fj) , and lGarciaetai] 42001b . The accretion model shows im- 
proved, though not acceptable, O-C agreement. 

ing an overshoot parameter of the order of 0.15 times 
the pressure scale height and increasing the opacity 
by ~7% (within the uncertainty limits of the abun- 
dances), they were able to reproduce the seismically 
inferred CZ base and Y cz . However, large sound- 
speed discrepancies remain in the radiative region of 
their model. 

To explore the possibility of convective overshoot, 
we evolve models with AGS05 abundances but extend 
the CZ that follows the adiabatic gradient to a depth 
that optimizes agreement with the sound speed inver- 
sions. We hoped that a deeper CZ would also inhibit 
diffusion and keep the CZ Y abundance higher. 

5.2. Results 

Figure QT| shows the relative sound-speed differ- 
ences for models with the GN93 mixture, the AGS05 
mixture, and the first convective overshooting model. 
For the first overshoot model, the CZ depth is 0.704 
R , deeper than inferred seismically. The sound speed 
agreement is improved only within the CZ but not 
much below it. The deeper CZ does not inhibit Y diffu- 
sion as we had hoped. In the second model, we extend 
the CZ even deeper, to 0.64 R . Figure fT2l shows the 
relative sound-speed difference for the second convec- 
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Fig. 11. — Relative difference between inferred and calculated 
sound speeds for models with the GN93 and AGS05 abundances and 
the first convective overshoot model that extends the adiabatically 
stratified layer to 0.704 Rq. 



Fig. 12. — Relative difference between inferred and calculated 
sound speeds for models with the GN93 and AGS05 abundances and 
the second convective overshoot model that extends the adiabatically 
stratified layer to 0.64 Rq. 



five overshooting model and standard models with the 
GN93 and AGS05 mixtures. The sound speed gradi- 
ent at the base of this adiabatically stratified CZ clearly 
does not agree with the seismically inferred one, and 
sound speed agreement in the central 0.2 R is much 
worse than in any of the other models. We therefore 
omit the second overshoot model from further analy- 
sis. Figure Q~3] shows the observed minus calculated 
frequency for the first overshoot model compared to 
the GN93 and AGS05 models. The AGS05 model and 
the overshoot model where the core Z is low do not 
agree as well with the data as the GN93 model does. 
It appears that overshooting alone is not a solution to 
this problem. 

6. Summary of mitigation attempts 

Table [2] summarizes the CZ Y, CZ base radius, and 
photosphere Z/X for the models examined here. Figure 
[T4l shows the relative sound-speed differences of our 
models with the GN93 mixture, the AGS05 mixture, 
and models with mass-loss, accretion, and convective 
overshoot. The observed minus calculated frequencies 
of these models are shown in Figure Q3] As seen in 
Tableland Figures Qj] and [H the AGS05 model and 
overshoot model do not agree well with the data. The 
GN93 model, the 1.3 M mass loss model, and the ac- 
cretion model show better agreement, though no model 
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Fig. 13. — Observed minus calculated frequency versus cal- 
culated frequency for models with the GN93 and AGS05 abun- 
dances and for the overshoot model for degree 1=0, 2, 10, and 20 
modes. The calculated frequencies were computed using the lPesneTI 
( 1990) non-adiabatic stellar pulsation code, and the d ata are from 
IChaplin et at] Jl998l) . ISchou & Tomczvld 119961) . and iGarcia et alj 
42001b . 



matches the data perfectly. 

Figure[l6]shows the small frequency separation dif- 
ferences of the Guzik et al. models minus the solar- 
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Table 2: Initial mass and surface abundances, mixing length parameter, and final abundances and CZ base for our solar 
models. 



Model 


GN93 


AGS05 


ML 1 


ML 2 


Accretion 


Overshoot 


C 1 


C2 


C3 


M o /M 


1.00 


1.00 


1.30 


1.15 


0.98 


1.00 


1.00 


1.00 


1.10 


Y 


0.26930 


0.25700 


0.24659 


0.25279 


0.26927 


0.25699 


0.26370 


0.27780 


0.27530 


Z 


0.01970 


0.01350 


0.01351 


0.01351 


0.01973 


0.01351 


0.01740 


0.01700 


0.01700 


a 


2.0379 


2.0004 


2.0571 


2.0104 


1.8958 


1.9962 


1.9918 


2.0635 


2.0652 


Z/X 


0.0240 


0.0163 


0.0178 


0.0171 


0.0162 


0.0164 


0.0209 


0.0208 


0.0219 


Ycz 


0.2412 


0.2273 


0.2388 


0.2349 


0.2402 


0.2292 


0.2349 


0.2473 


0.2551 


Rczb /R© 


0.7125 


0.7294 


0.7217 


0.7264 


0.7241 


0.7038 


0.7186 


0.7190 


0.7181 



"Seismically inferred values from lBasu & Antial i2004): Y cz = 0.2485 ± 0.0035, R C zb /Rd = 0.713 ± 0.001. 

''Models ML 1 and ML 2 are the AGS05 models including mass loss. Models CI, C2, and C3 are the C0 5 BOLD models with GN93 opacities, 
AGS05 opacities, and AGS05 opacities including mass loss. 
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Fig. 14. — Relative difference between inferred and calculated 
sound speeds for models with the GN93 and AGS05 abundances 
and models with mass-loss, accretion, and convective overshoot. 



cycle corre cted frequency diffe rences from the Bi- 
SON group dChaplin et al.ll2007l) for ^=0 and 2 modes, 
which are sensitive to the structure of the core. This 
plot illustrates that including low-Z accretion in the 
model retains the core structure of the GN93 model. 
The overshoot model and AGS05 model do not agree 
as well with the data as the models with higher core Z. 
The mass losing model with the best sound speed and 
observed minus calculated frequency agreement (Mo = 
1.3 M Q ) shows worse agreement than any of the other 
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Fig. 15. — Observed minus calculated frequency versus calcu- 
lated frequency for degree (=0, 2, 10, and 20 modes in GN93 and 
AGS05 models and in models with mass-loss, accretion, and con- 
vective overshoot. The calculated frequencies were computed us- 
ing the Pesnell ( 1990) non- adiaba t ic stellar pulsation code. The 
data are fromlChaplin et al] <2007h ,l Schou &Tomczvl<M 19961) , and 
iGarciaetallfeOOll) . 



models. However, mass loss does change the value of 
the small frequency separation in the right direction to 
correct for the discrepancy seen with the new abun- 
dances. Perhaps the over-compensation indicates that 
this model has too much mass loss; a model with a 
smaller initial mass would reduce the disagreement, as 
seen with the M = 1.15 M G model. 
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Fig. 16. — Difference between calculated and observed small sep- 
arations for 1=0 and 2 modes for the Guzik et al. models. The data 
are from Chaplin et ai] j2007l) . 

Mass-losing models can improve seismic agree- 
ment for the new abundances, but they do not fully 
restore agreement. In addition, the destruction of too 
much Li and the production of too much surface 3 He 
make the two models considered here unlikely. A 
smaller amount of mass loss that leads to destruction 
of some, but not all, of the initial Li could provide 
a plausible partial mitigation of the solar abundance 
problem. 

The accretion model allows for a solar interior that 
is similar to models developed with the higher abun- 
dances and therefore agrees nicely with seismic in- 
ferences in the central 0.5 R of the sun. The model 
still shows poor agreement near the CZ base. The very 
steep Z abundance grad ient develo ped at the CZ base 
seen in Figure l7ksee iGuzikl 12006) might have a de- 
tectab le signature in the seismic frequencies. iBasu 
d 19971) finds that inversions appear to rule out such 
steep composition gradients at the CZ base. 

The convective overshoot model with a CZ depth 
of 0.704 R Q improves sound speed agreement slightly 
within the CZ but not much below it. Extending the 
CZ depth to 0.64 R results in even worse agreement. 
In addition, Y diffusion is not inhibited by the deeper 
CZ, as we had hoped. 



7. CCFBOLD abundances 

At the suggestion of both the referee and a col- 
league R Bonifacio, we also include here a prelimi- 
nary exploration of solar models using the CO 5 BOLD 
abundances with a Z/X of 0.0209 and Z of 0.0154, 
intermediate between the AGS05 and AGSS09 abun- 
dances, but lower than the GN93 abundances. 

Because the abundances of only 12 elements have 
been re-evaluated at this time by the CO s BOLD group 
dLudwig et al. I l2009h . it is a little premature to create 
new opacity tables for the CO s BOLD mixture, which 
can, in principle be done using the Lawrence Liver- 
more OPAL web request at http://opalopacity.llnl.gov 
In addition, we do not have low-temperature opacities 
for a mixture representative of the CO 5 BOLD abun- 
dances available to us. Therefore, here we decided 
to calibrate standard models to the Z/X of CO 5 BOLD 
using opacity tables based on the GN93 and AGS05 
mixtures. We observe that the O/Fe (oxygen to iron) 
mass fraction ratio of CO s BOLD is 4.98, intermedi- 
ate between the O/Fe mass ratio of 6.56 for the GN93 
mixture and 4.65 for the AGS05 mixture, so our two 
standard models should bracket results that use the 
CO s BOLD mixture in the opacity tables. We did up- 
date abundances to the CO s BOLD values for our in- 
line equation of state calculation and for tracking the 
diffusion of the major elements. 

The initial Y, Z, and mixing length to pressure- 
scale-height ratios needed to calibrate to the CO s BOLD 
Z/X using either opacity set are listed in Table [2] Fig- 
ures [17] [TU and [19] show the results for sound-speed 
differences, small separations between calculated f-Q 
and 1-2 modes, and observed minus calculated fre- 
quencies for ( = 0, 2, 10, and 20. The sound speed dis- 
crepancy is reduced to only about 0.6% at the convec- 
tion zone base for the CO s BOLD abundances, com- 
pared to 1.4% for the AGS05 abundances and 0.4% 
for the GN93 abundances. The results for either opac- 
ity set are identical above 0.6 R Q , but differ below this 
where oxygen is the main opacity contributor, and in 
the core where iron is the main opacity contributor, as 
expected. The model using the AGS05 opacity mixture 
requires a higher Y abundance to compensate for the 
relatively higher Fe opacity contribution in the core. 
The small separations and observed minus calculated 
frequencies are slightly higher than found for a model 
calibrated to the GN93 abundances on average, but not 
as high as for a model calibrated to the AGS05 Z/X. 

As also surmised by the referee, since the CO s BOLD 
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sound speed differences are closer to observed, im- 
provement could be obtained with a smaller amount 
of mass loss. Here we have calculated an addi- 
tional mass-loss model with initial mass 1.1 M and 
an initial mass-loss rate of 2.25 x 1CT 10 M /yr, ex- 
ponentially decaying with e-folding time 0.45 Gyr. 
This model was calibrated to the CO 5 BOLD Z/X us- 
ing the AGS05 opacities that have O/Fe abundance 
close r to that of the CO 5 BOLD O/Fe . Previous work 
(e.g. ISwenson & Faulknerill992l: iGuzik & Cox|[l995 



Sackmann & Boothrovdl 2003b indiciated that an ini- 
tial mass of 1 . 1 M or less and a relatively short mass 
loss phase (less than 0.2-0.5 Gyr) could deplete the 
lithium to the present-day observed values from initial 
solar-system abundance without completely destroy- 
ing the lithium or building up too much 3 He. The 
sound speed agreement (Figure ITTb shows consider- 
able improvement with this smaller amount of mass 
loss; however, the agreement for the solar core is not 
as good as for the non-mass losing models, as can be 
seen more clearly in the small separations (Figure [18]). 
This mass-losing model restores the level of agreement 
attained with the GN93 abundances for the observed 
minus calculated frequencies (Figure [T9l. 

Figure|20]shows the luminosity versus time for this 
model, as well as for the standard solar models evolved 
with GN93 and AGS05 abundances. Figure |2TI shows 
the effective temperature experienced by the surface 
layer throughout the evolution of the mass-loss model 
compared to that experienced by the standard models. 
During the mass-losing phase, the high temperatures 
that depleted the lithium in the current surface layer 
were experienced when the material that is at the sur- 
face of the now 1 M sun was deeper, before the previ- 
ous surface layers were lost. After the mass-loss phase, 
the temperature that affects the surface layer is that of 
the CZ base, since the material currently at the surface 
is continually mixed through the CZ. We see that in- 
cluding mass loss in the CO s BOLD model exposes the 
surface layers to high enough temperatures to deplete 
Li early in the evolution (2.8 million K is the temper- 
ature required for relatively rapid Li destruction). It is 
not clear that the advantages of mass loss (e.g. Li de- 
pletion and better sound speed agreement in the outer 
80% of the solar radius) can be retained while at the 
same time not creating a discrepancy in the inner 20%. 
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Fig. 18. — Difference between calculated and observed small sep- 
arations for £=0 and 2 modes for the Guzik et al. models using the 
CO 5 BOLD abundances with either the GN93 or AGS05 opacities 
and with the AGS 05 opacities and mass loss. The data are from 
IChaplin et aT] <2007l) . 



8. Conclusions and future work 



In spite of the seismic evidence in favor of the old 
abundances, the new abundances cannot be easily dis- 
missed. The improvements in the physics of the atmo- 
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spheric models used to determine the abundances, the 
success achieved in line-profile matching, and the self- 
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ties. For the mass-losing phase, the lithium-destroying temperatures 
are attained because the layer that is now at the surface once resided 
deeper inside the sun. In the post-mass-loss phase, the relevant tem- 
peratures are attained by envelope convection which mixes surface 
layers downward, exposing the surface material to the temperature at 
the CZ base. 2.8 million K is the temperature required for relatively 
rapid Li destruction. 



consistency of the abundance determinations provide 
great credibility to the new lower abundances. How- 
ever, solar models developed with the new abundances 
remain discrepant with helioseismic constraints, even 
with a variety of (often unjustified) changes to the in- 
put physics. Adjustments to the evolution of solar 
models, such as the early mass loss and low-Z accre- 
tion discussed here, show some promise but do not 
fully restore agreement. Any single adjustment to so- 
lar models does not fully resolve the problem. Combi- 
nations of changes might provide better agreement but 
seem contrived. A resolution to the solar abundance 
problem (or the solar model problem) remains elusive. 

In the future, a more comprehensive exploration of 
parameter space, including opacity variations, differ- 
ent mass-loss or accretion prescriptions, diffusion, and 
perhaps even combinations of these effects could be 
useful. In particular, models with AGS05 abundances 
and a smaller amount of mass loss than explored here 
may provide a way to retain the core structure with- 
out completely destroying Li or creating too much 
3 He build-up. Further examination of the C0 5 BOLD 
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models with revised abundances for every element and 
opacity tables (including low-T opacities) adjusted for 
a new mixture would also be enlightening. 
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